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SUMMARY 

The first part of this technical note covers the study of 
guidance of HF radiowaves along field-aligned paths using the top- 
side sounder data from the Alouette 2 satellite. More than 100,000 
Alouette 2 ionograms have been used in this analysis and an un- 
usually high percentage of occurrence of conjugate echoes has been 
observed in the data of some equatorial stations. The maximum per- 
centage of occurrence of conjugate traces around local midnight 
when the satellite is near the apogee (3000 km) throws some light 
on the diffusion mechanism responsible for the formation of field- 
aligned irregularities. The Singapore data (Geomagnetic longitude 
180° ±30°) has shown higher occurrence of ducted echo-traces (>10 
percent) than the data from all other stations. This is rather a 
large percentage. The longitudinal preference may be conducted 
with the larger magnetic field strengths in that longitudinal range 
Two interesting types of ionograms containing conjugate echo traces 
are presented. They are called the "double hook" and the "triple 
hook" types. A self-consistent explanation based on the assumption 
of multiple reflections between conjugate points of the field line 
passing through the satellite is presented. 

The second part of the note covers the results of the digital 
ray-tracing study of the guidance of HF radiowaves along field- 
aligned ionization irregularities in the magnetosphere of the Earth 
Various computer results are shown to illustrate the mechanism of 
guidance of radiowaves. Time-delay and trapping characteristics 
of the two types of conjugate echoes that are frequently observed 
in the Alouette 2 topside sounder records are simulated using the 
ray-tracing methods . Some of the optimum conditions for guidance 
of electromagnetic waves along a magnetic field line are discussed. 
The results reported here on the HF ducting phenomenon is directed 
primarily towards determining whether such a mechanism may be 
utilized to study the magnetospheric phenomena and whether it can 
be useful as a reliable communication link. 


*The information in Part I of this report is also contained in a paper entitled '’Field-Aligned Propagation HF 
Radiowaves Using Alouette-2 Topside Sounder Data,” by J. Ramasastry, E. J. Walsh, and J. R. Herman; 
and the material in Part II is contained in a paper entitled "Field-Aligned Propagation HF Radiowaves Using 
Digital R.ay-Tracing Techniques,” by J. Ramasastry, E. T. Walsh, and J. R. Herman. Both of the papers 
were presented at the AGARD Meeting, August 19-23, 1968, in Oslo, Norway. 



PART I 


CONJUGATE ECHOES IN ALOUETTE-2 TOPSIDE SOUNDER DATA 


Introduction 

High frequency conjugate echoes were first observed on 
the Alouette 1 sweep frequency topside sounder data (ref. 1) . 
Conjugate signals have also been observed on the Explorer 20 
fixed-frequency topside sounder ionograms (ref. 2). Some 
occurrence statistics concerning the Alouette 2 ionograms have 
also been published (ref. 3). However, the study of HF conjugate 
echoes on Alouette 2 ionograms has offered challenging problems 
and some exciting results. Two types of ionograms containing 
symmetric conjugate echo traces which appear frequently are pre- 
sented and explained in terms of multiple reflections between 
the conjugate points of a field line passing through the satel- 
lite. Such symmetric traces have not been previously reported 
in the literature. 

The Alouette 2 topside sounder is better suited to the 
study of conjugate echoes than that of Alouette 1 because of the 
lower frequency limit of the sounder. The transmitter and the 
receiver are tuned to the same frequency which is swept in time. 
The transmitter is pulsed once in 32 msec and the instrumentation 
provides a new zero time-reference each 32 msec, thus pulses 
returning with delays of 37, 69, or 101 msec will all be recorded 
as traces of 5 msec group-delay on the same ionogram frame. 

More than 100,000 ionograms from eight stations (Quito, Fort 
Myers, Santiago, Johannesburg, Kano, Singapore, Winkfield, and 
Orroral Valley) have been examined by the authors and some of the 
observations of the HF ducting phenomena in the magnetosphere of 
the Earth have been very revealing. Various statistical results 
are presented and the occurrence characteristics of the conjugate 
echoes are discussed. 


Method of Analysis 

All the available ionograms from the eight stations were 
analyzed to identify the ionograms with conjugate traces. 

Muldrew (ref. 3) classified all the available ionograms into 
five groups depending on the value of f X s» the frequency of the 
extraordinary wave which has zero refractive index at the satel- 
lite altitude. The reason given by the author for such a 
classification is that it gives a rough indication of the elec- 
tron density at the satellite altitude. While it is a convenient 
method of classif ication , it does not provide any new physical 
insight into the guidance of HF radiowaves. Consequently, the 
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authors have not adopted any method of classification of frames 
with conjugate echoes. The statistical data in the present 
analysis are restricted to those ionograms which exhibit con- 
jugate echoes. Hence the percentage occurrence should be taken 
to mean that so many out of 100 ionograms with conjugate ducts 
occur at the indicated value of the parameter (altitude, local- 
time, or L-value) . 

Since the argument that the perigee rotates only 1.89 deg/ 
day, which is almost the same as the precession of the orbital 
plane with respect to the sun, the local time at perigee (or 
apogee) changes less than 3 hr/yr. For the period of the data 
used in our analysis, the perigee occurs near local noon and the 
apogee near local midnight. Hence the satellite intersects the 
smaller L shells near noon and the larger L shells near midnight. 
As we shall see later, the percentage occurrence is rather large 
around local midnight when the satellite is near the apogee. 


Conjugate Echoes as a Function of L Value 

The percentage occurrence of conjugate echoes between 
21 hours and 04 hours local time is plotted as a function of L 
in Figure 1. This period of local time is chosen to afford 
comparison with the results of Muldrew (ref. 3). The results 
presented here are consistent with Muldrew' s argument that the 
occurrence percentage is high and almost constant during this 
time period. The results for L values larger than 3.5 are not 
considered representative and hence may be deleted from discus- 
sion. Very few conjugate ducts are observed at these high L 
values because the receiver is no longer tuned to the frequency 
of echoes that arrive after the long delay required for propaga- 
tion along a duct with an L value about 4 or greater. Such a 
conclusion agrees with our analysis of "triple hook" and "double 
hook" traces where the traces corresponding to multiple reflec- 
tions between the conjugate point start fading out as the group 
delay exceeds 200 msec. The authors' results also show an 
apparent drop in occurrence for L less than about 1.35. It can 
also be agreed that such a rapid drop may not be real since the 
data at such low L values are not representative because of the 
characteristics of the satellite orbit, i.e., the satellite 
intersects lower L shells during local noon where the percentage 
occurrence of conjugate echoes is almost zero (ref. 3). One may 
say in this context that the Alouette 2 satellite was not designed 
originally to study conjugate ducting. Future experiments 
designed to study conjugate ducting should have the capability of 
receiving signals for longer periods of time at any one frequency; 
in addition, the orbital characteristics of the satellite should 
be such that the data are representative of guidance in field 
lines of both low- and high-L values. 
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Figure 1.- Percentage occurrence of conjugate ducts as a function of L value 



Conjugate Echoes as a Function of Local Time 
(1-hour averages) 

The percentage occurrence of conjugate echoes as a function 
of local time is shown in Figure 2. The data were averaged in 
1-hour intervals of local time. This averaging time interval is 
selected in view of the fact that the local time drift along the 
propagation path may be about an hour. It can be seen that the 
percentage occurrence of conjugate echoes is almost zero in the 
time interval 0600 to 1800 hours local time. There is a pro- 
nounced broad maximum between 2100 and 0400 hours local time. 

The post-midnight maximum in duct occurrence when the satellite 
altitude is near apogee brings forth the discrete characteristics 
of the field-aligned ionization irregularities. As may be noted, 
the ambient density during night time is considerably lower and 
any diffusion mechanism capable of channelling the particles 
along the field line creates larger gradients of electron density 
along the field line during the night time. Hence the increased 
efficiency in the guidance of radiowaves during night time can be 
anticipated. 


Conjugate Echoes as a Function of Satellite Altitude 


The percentage occurrence of conjugate echoes as a function 
of satellite altitude is shown in Figure 3. As mentioned before, 
the satellite is around the midnight meridian when it approaches 
the apogee of its orbit. Hence the high altitude maximum is con- 
sistent with the previous conclusions about midnight maximum. 
Conjugate echoes are rarely seen at altitudes less than 1200 km 
in the Alouette 2 data. 


Examination of Conjugate Traces 


In the Alouette 2 ionograms, the conjugate echoes generally 
appear as beautifully symmetric patterns of high intensity. When 
they appear, the normal traces in the ionogram are faded out be- 
cause of the automatic gain control in the receiver. Figure 4 
shows a series of evenly spaced "double hooks" separated by hori- 
zontal traces, while Figure 5 shows a pattern of symmetric 
"triple hooks." The symmetry of these patterns is in sharp con- 
trast to the types of conjugate traces reported by Muldrew 
(ref. 1) for Alouette 1 data. They are also in sharp contrast 
to the conjugate traces reported by Muldrew (ref. 3) for Alouette 
2 data. The L values of the magnetic field lines passing through 
the satellite for these two examples (L = 1.53, 1.63, respec- 
tively) are larger than the L values for the conjugate echo traces 
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LOCAL TIME (HOURS) 

Figure 2.- Percentage occurrence of conjugate ducts as a function 
of Local Time (1-hour average) 



ALTITUDE (KMS) 

Figure 3.- Percentage occurrence of conjugate ducts as a function 
of altitude (100-km average) 



0.9 123 IS 20 NORTH | 

FREQUENCY (UH«) _»„>>*>• ' »• 


Figure 4. - The "double hook" pattern recorded 
at Quito, Ecuador, on Aug. 18, 1SF66, at 
08:00:05 GMT. (Local Time = 02:07:46, Geoma 
Lat. = 12.56°N, Geomag. Long. = 18:22°W, 

L = 



FRKHJEWcr ( MM* ) 


Figure 5. - The "triple hook" pattern recorded 
at Santiago, Chile, on April 22, 1966, at 
09:01:48 GMT. (Local Time = 04:38:41, Geomag. 
Lat. = 1 9 . 53°S , Geomag. Long. = 2.54°E, L = 1.63) 
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considered by Muldrew (ref. 1) for the Alouette 1 data. Hence, 
the major portions of the paths are in the regions of low elec- 
tron densities for the cases considered here. This means that 
the group delay for a round trip path (N + S) should be rela- 
tively independent of frequency and approximately determined by 
the geometrical length of the field line between the two con- 
jugate reflection points. N and S correspond to the north and 
south directions taken by the signal. Only the extraordinary 
wave traces of conjugate echoes are observed. Hence all the 
round trip traces (N + S, 2N + 2S, 3N + 3S,...) should have a 
tendency to appear as horizontal lines in the ionograms. The 
symmetry of the "double hook" and "triple hook" patterns follows 
directly from this characteristic. The far end echo from the 
satellite always heads along the field line away from the Earth 
towards the region of decreasing electron density, except at the 
equator where the path-lengths and the characteristics are the 
same. Hence, reflections were not obtained from near the satel- 
lite as for the case of the near-end echo which follows the field 
line towards the Earth and towards the region of increasing elec- 
tron density. This results in reflections from near the satellite. 
The subsequent symmetry of the N and S traces follows from the 
magneto-ionic properties of the ray path. 

The ionogram in Figure 4 is a typical example of the fre- 
quently observed "double hook" patterns. It was recorded at 
Quito, Ecuador. Also shown in the figure are the proper rela- 
tive positions of the various traces according to their actual 
time delays. The satellite position (GM Lat = 12.56°) is such 
that the group-delays of the near-end and far-end paths along the 
field line differ approximately by 32 msec and the North and South 
traces seem to intersect in real time. The geometry of the path 
length of the field line requires the round trip trace (N + S) to 
be positioned twice as far from the beginning of its own time 
frame as the apparent intersection point of the N and S traces. 

The (2N + 2S) trace should be twice as far from the beginning 
of its frame as the (N + S) trace is from the beginning of its 
own time frame. Similar geometric considerations can account 
for the uniform spacing of the traces in the ionogram. The 
(2N + S) and (3N + 2S) traces are shaped like the N trace, and 
the (N +2S) and (2N + 3S) traces are shaped like the S trace. 

This is because the delay times for the round trip traces (N + S) , 
(2N + 2S) , (3N + 3S),... are virtually independent of frequency. 

The direct consequence of such a situation is that all the 
"double hook" traces do have the same shape. The (3N + 3S) trace 
is extremely faint but is discernible on the microfilm. 

Figure 5 shows the "triple hook" ionogram and its various 
characteristics. The reason that all the triple hooks in the 
ionogram of Figure 5 have the same shape is once again caused by 
the frequency-independent group-delay characteristics of the 
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round trip traces. The figure also shows the various traces 
according to their respective time delays. The near-end trace 
is easily recognizable in both Figures 4 and 5 because it is the 
most intense trace which starts from zero group delay. This 
characteristic is caused by reflections from near the satellite 
due to the increasing electron density in that direction. In 
both the ionograms, the far end trace was initially selected on 
the basis of its intensity and the mirror-image-like character- 
istic with respect to the near-end trace. Also, the L value of 
the field line corresponding to that frame gave the clue as to the 
approximate group-delay of the far-end echo. Digital ray tracing 
analyses have, since then, borne out these assumptions. 

Figure 5 exhibits triple hooks instead of double hooks 
because the satellite is displaced from the magnetic equator 
(GM Lat = -19.53°) so that the group-delays of the near-end and 
far-end echo traces differ by more than 32 msec. The apexes of 
all the triple hooks occur at the frequency where the time delay 
of the near end (S) echo is exactly one frame time (32 msec) . 

Then the N, N + S, and N + 2S echoes all have the same positions 
in their relative time frames at that frequency and will seem to 
intersect in the ionogram where all the time frames are super- 
imposed . 

There is one significant condition which determines whether 
the conjugate echoes form a triple hook pattern or a double hook 
pattern. Double hooks are observed when the group-delays of the 
near-end and far-end echoes differ by approximately 32 msec. 

Triple hooks are observed when the group-delays differ by more 
than 32 msec. 

The traces representing the round trip paths are sometimes 
inclined to the horizontal. The slopes of (2N + 2S) , (3N + 3S) , 

(4N + 4S) , and so forth, are exactly twice, thrice, four times, 
and so forth, the slope of the trace (N + S) . This can be easily 
verified by adding up the various group delays. Such inclined 
round-trip traces appear when the near-end and far-end distances 
differ by a small amount. 

One of the significant features observed in Figure 5 is 
that the (4N + 4S) trace had a delay time of 465 msec. But with 
the receiver bandwidth of 30 kHz and the sounder frequency-sweep 
rate of 0.15 MHz/sec, the receiver is essentially detuned after 
200 msec. This may mean that the faintness of the (4N + 4S) 
trace is not caused by the leakage and attenuation of the signal 
within the "magnetospheric wave-guide." On the contrary, the 
signal after transversing four complete round trips from one 
hemisphere to the other, covering a total path-lenght of about 
60,000 km, must have been quite strong to have shown up at all 
in the detuned receiver. 
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The faint (3N + 3S) trace in the "double hook" case 
(Figure 4) has a time delay of 316 msec and probably would have 
been more intense if the receiver had been properly tuned to 
receive it. 

Figure 6 shows a conjugate echo frame corresponding to a 
high L value of 1.76. The satellite was at -19.05° (mag. lat.) 
and the total field-aligned path was 17,463 km with the near-end 
and far-end distances being 6593 and 10,870 km, respectively. 

Figure 7 shows a series of seven ionograms recorded at 
Quito on June 12, 1966, from 03:23:45 to 03:26:59 UT. As may 
be seen the triple hooks prominently appear and disappear. Such 
a sequence is representative of an event when the satellite 
passes through a cluster of ionization irregularities with strong 
coupling of signals into ducts of various closely spaced L shells. 

One may also emphasize in this context that the ionization 
irregularities transverse to the field lines responsible for such 
a wave guidance may be either enhanced or depleted columns of 
ionization. The nature of the irregularities is not discussed 
in the present note. 


Discussion and Conclusion 

One must bear in mind that the Alouette 2 satellite topside 
sounder experiment was not originally designed for detecting such 
long time-delay echoes. The frequency sweep rates and bandwidths 
used are essentially sufficient to receive the signals normally 
incident on the topside ionosphere. The Explorer 20 satellite 
(ref. 2) had a fixed frequency topside sounder transmitting at 
seven different frequencies in sequence, each time the cor- 
responding receiver being tuned for 15 msec. The next time the 
receiver could receive the signal at that frequency was 105 msec 
later. Since the satellite altitude was only 1000 km (circular 
orbit) like that of Alouette 1, the experiment was once again 
not suitable for the study of conjugate ducting. A suitably 
designed satellite propagation experiment, with its receiver 
being capable of receiving signals of larger time delays, may 
unravel all the mysteries about magnetospheric wave-guidance. 

The authors also propose an experiment in which a trans- 
mitter in a low orbiting satellite (preferably eccentric polar 
orbit) transmits at HF to a high-altitude satellite (altitude of 
apogee ~6000 km) whose receiver should have the capability of 
receiving the signals for longer periods of time. The now 
orbiting Radio Astronomy Satellite might have the desired re- 
ceiving characteristics in the frequency range 1-10 MHz. Many 
monostatic and bistatic radar experiments of the past aimed at 
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Figure 6. - Ionogram showing conjugate echoes for a high L value (L = 1.76) 
recorded at Johannesburg, South Africa, on Aug. 12, 1966, at 02:02:54 GMT 
(Local Time = 02:37:05, Geomag. Lat. = -19.05°S, Geomag. Long. = 74.53°E) 
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Figure 7. - A series of seven ionograms recorded 
at Quito, Ecuador, on June 12, 1966, from 
03:23:45 to 03:26:59 GMT. 
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verifying the existence of ducts and shells in the magnetosphere 
and the consequent guidance of electromagnetic waves have been 
inconclusive . The stability and frequency of occurrence of the 
conjugate echoes combined with the relatively loss-free multiple 
reflections between the conjugate points, of the signals observed 
in the Alouette 2 satellite data have created challenging oppor- 
tunities. Thus a well designed experiment may provide conclusive 
answers to the question of whether or not this unique magneto- 
spheric wave guide could be used as a reliable HF communication 
link between conjugate points or between two satellites over 
large transequatorial distances along a field-aligned path. Even 
if this propagation mode cannot be used as a communication link, 
its use to study the magnetospheric phenomena cannot be dis- 
counted . 
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PART II 


DIGITAL RAY-TRACING STUDIES 


Introduction 

Some of the echo traces appearing at virtual ranges greater 
than those of the normal vertical incidence echo traces on the 
topside sounder ionograms have been explained in terms of field- 
aligned ionization irregularities (ref. 1). The irregularities 
are assumed to have thicknesses greater than the radio wavelength 
and to act as waveguides or "ducts" to trap HF energy and produce 
the long-range echo traces. A complete review of the past 
theoretical and experimental work directed towards testing such 
a guidance mechanism has been published (ref. 4). 

The guiding of rays along field lines requires irregular- 
ities with a certain minimum transverse ionization gradient. 
Propagation between magnetic conjugate points or "conjugate 
ducting" occurs when the transverse ionization gradient exceeds 
the minimum required for guidance at the apex of the magnetic 
line of force. The current interest in the HF ducting phenomenon 
is directed primarily towards determining whether such a mecha- 
nism may be utilized to study the megnetospheric phenomena and 
whether it can be useful as a reliable communication link. Here, 
we also examine the conditions under which the HF radiowaves may 
be guided along field-aligned ionization irregularities. 


Ray-Tracing Method 

The ray tracing method is based on Hamilton's system of 
equations (ref. 5) as derived in spherical polar coordinate 
system by Haselgrove (ref. 6) and extended by Grossi and 
Langworthy (ref. 7) for the investigation of HF and VHF iono- 
spheric progagation. The reader is referred to reference 7 for 
a description of the basic ray-tracing program. Because of the 
high accuracy obtainable, the use of a high-speed digital com- 
puter is preferred to an analog machine for the integration of 
the ray equations. The Adams-Moulton variable step-size inte- 
gration method (ref. 8) is adopted because of its greater flexi- 
bility due to an automatic error control which halves or doubles 
the step size to ensure the desired accuracy along all parts of 
the ray path. 


Magnetospheric Models 

The electromagnetic properties of the magnetosphere are 
characterized by (1) the electron density distribution, (2) the 
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magnetic field distribution, and (3) the electron collision- 
frequency distribution. The mathematical models for the above 
three items are implemented in the form of subroutines and can 
therefore be modified as and when necessary. 


Magnetic Field Model 

To a first approximation, a dipole model is used for the 
magnetic field of the Earth. The magnetic field equation which 
defines the gyrof requency f H is given by: 


f H (r,0) 



3cos^ 0 


( 1 ) 


where a is the radius of the Earth, and r and 0 are the geo- 
centric distance and colatitude of any point on the field line. 


'11 


B x 10 
2Trm o 


-6 * 


0.9 


( 2 ) 


where B 0 is the magnetic field on the surface of the Earth at the 
equator, and e and m are the charge and mass of an electron. A 
value of 0.3142 gauss is used for B 0 . 

For ray tracing along high L-value field lines, the 63- 
term spherical harmonic expansion of Cains, Daniels, Hendriks, 
and Jensen (ref. 9) is used. 


Collision Frequency Model 

The collision frequency model is based on results published 
by F.S. Johnson in the Satellite Environment Handboo k (ref. 10). 


Electron Density Model 

The model for the electron density distribution incorporates 
those features that have a bearing on the high-frequency ducting 
problem. Thus, an ionization irregularity normal to the field 
line is superimposed on the normal radial distribution of elec- 
tron density. In this model, the electron density is given by: 

N (r,0 , <J>) = N 1 (r,0)*N 2 (r,0,<|)) (3) 
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where 0 is the colatitude and 4> is the magnetic longitude. 

N^(r,0) represents the background density in the magnetic meri- 
dional plane and N2(r,0,4>) represents the ionization in the 
field-aligned irregularities in both the meridional and azimuthal 
planes . 

The background ionization is assumed to have the form: 


N-^r,©) = N p + N x . 


(4) 


Here Np and Njj refer to the electron densities in the ionospheric 
F-region and the exosphere, respectively. 


N 


F 


N 

max 


*exp 




(5) 


where N is the peak electron density and 

max J 


W 


h - h 


max 


H„ 


(6) 


in which h max is the altitude of the peak electron density, N max 
and Hp is the scale height. In our calculations, the maximum 
altitude h max is assumed to be 350 km and Hp = 50 km. 

The Angerami and Thomas diffusive equilibrium model (ref. 
11) is used for the exospheric electron density distribution 
above 1000 km and below the plasma-pause which exists at a dis- 
tance of from 3 to 4.5 Earth radii. The centrifugal force 
caused by the rotation of the Earth is neglected, while applying 
the Angerami and Thomas model. With this model, the electron 
density in the exosphere is given by: 


N x = N r * 


t 1 ’ 


659 x 


1O" 3 0 


+ 1.662 x 10 


*] 


(7) 


where 0 is the colatitude in degrees. The terms within the 
brackets represent the latitudinal variation of exospheric elec- 
tron density as derived by Thomas and Sadar (ref. 12) using 
Alouette 1 satellite data. 
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( 8 ) 


N (Z) 
r 



where Z Y is the geopotential altitude given by 


Z X - ^ (r - r o> (9 > 


where r 0 is the reference radial distance and is equal to 
6870 km in the model. The symbol i refers to each ion present 
in the exosphere. Hi and 6i refer to the scale height and the 
fractional density of the i^h ion at the reference level r Q . 
The scale height is given by 


kT . 

H = 1 - 
i M ± g(r 0 ) 


( 10 ) 


where k is the Boltzmann Constant, Ti and Mi are the temperature 
and mass, respectively, of the i th ion, and g(r Q ) is the acceler- 
ation due to gravity at the reference level r 0 . A three-ion 
gas model consisting of oxygen (0 + ) , helium (He + ) , and Hydrogen 
(H + ) ions is assumed. The various parameters in the density 
equation have the following values: 


T = 1000°K for both electrons and ions 


0.2 x 10 


-T/1000 


He 




= 1.96 x 10 


-2 


Bh = 0-16 X^10- T / S0 ° _ 16 x 10 -3 


B 0 - t - °- 9788 

e 


(Ha) 

(lib) 

(He) 

(lid) 
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B e = 1 + 0.2 x 10 (- T /1000) + 0>16 x 1Q ( T/500) 


= 1.0216 


(lie) 


The technique suggested by Swayze (ref. 13) is used to 
join the exospheric profile smoothly with that of the F-layer 
between 350 and 1000 km. The expression for N x is modulated by 
a factor 

-p[-(^4o^) 2 ] 

where h is the altitude in kilometers. 

The electron density as a function of altitude at various 
geomagnetic latitudes is shown in Figure 8. 


Model for the Ionization Irregularities 

The following assumptions are made in the development of 
this model: 

1. The ionization in the irregularities is distributed 
in discrete columns or sheets separated in latitude 
and aligned with the magnetic field lines. 

2. The transverse scale of the field aligned irregu- 
larities follows the transverse (meridional) dimen- 
sions of the tube of magnetic flux. 

3. The peak enhancement varies inversely as the trans- 
verse dimensions of the tube of magnetic flux. 

4. The variation of the density within the irregularity 
follows a Gaussian distribution transverse to the 
magnetic field. The Gaussian distribution satisfies 
the equations for the diffusion of charged particles 
across a field line. 

Hence: 


N 2 (r , 0) 


1 + 



(13) 
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N is the peak ionization enhancement and H 0 is the scale size 
of the irregularity at the base of the field line. 


t _ Meridional wid th at (r, 0 ) 

t_ Meridional width at (r , 0 ) 

° o' o 

sin 3 0(4 - 3sin 2 Q )% 

~ 3 2 5T (14) 

sin -3 0 Q ( 4 - 3sin z 0) 2 


While specifying the initial coordinates of the starting 
point, the geographic latitude, A a , and the longitude, <j> a , are 
converted to geomagnetic coordinates A m and <j> m by using the 
following transformation equations: 


where 


and 


tan A m = -tan(A a - A Q )sin x sec (x + <J> ) 

tan = -cos A tan (x + <f> ) 
m m T o 

x = arctan-lcos (A - A ) cot d> X 
( a o T a( 

4> = +78.5° and A = +291.0° 

o o 


(15) 

(16) 


(17) 


are the geographic coordinates of the geomagnetic pole. 
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RESULTS 


Figures 9 through 13 show a few results of the ray tracing. 
Figures 9 through 12 simulate some of the conjugate echo paths 
observed in Alouette 2 ionograms . Parameters shown in the figures 
are defined as follows: 

R = the geocentric distance (km) of the initial 
° signal position, 

PKFRAC = the peak fractional ionization enhancement, 

LAMBDA = the colatitude of the field line, 

A = the initial angle the ray makes with the 
° local vertical, 

B = the initial angle the ray makes with the 
° South vector, 

PHI and THETA = the magnetic longitude and colatitude of the 

initial signal position. 

The frequency of the signal is 1.0 MHz for the cases considered 
in Figures 9 through 12. 

Each figure consists of two plots. The rectangular plot 
shows the distance of the ray from the field line versus the 
distance along the field line from the near end. The polar plot 
shows the actual ray path with reference to the surface of the 
Earth. The polar plots gives a clear indication of the positions 
of the conjugate reflection points as well as the L value of the 
field line guiding the rays. 

Figure 9 shows the L = 1.76 case using a peak electron 
density enhancement of 13.0 percent at the point the ray is 
launched. A peak enhancement of 13.0 percent at a geocentric 
distance of 10,022 km corresponds to a much smaller percentage 
at lower altitudes. The ray starts at about 12,500 km along the 
field line from the northern hemisphere. It gets reflected on 
its way towards the Earth and travels all the way back to the 
southern hemisphere where it is reflected again. The program 
was stopped deliberately when the ray was at 15,500 km along 
the field line. Otherwise, it could have bounced back and forth 
indefinitely. The bowing of the pattern at the equator is caused 
by larger scale sizes near the equator so that the ray rides 
further out from the field line because the enhancement structure 
is wider. 

The oscillatory nature of the ray path seen in the rec- 
tangular plot is caused by the ray being launched too far away 
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where the electron density gradient is insufficient to bend the 
ray back towards the direction parallel to the field line. The 
ray essentially tends to continue in a straight line, but appears 
to fall in towards the field line since the field line is curved. 
When the ray gets close enough, the gradient is large enough to 
bend the ray away from the field line and the apparent inward 
motion stops in the rectangular plot. It follows that similar 
oscillations also occur if the ray is launched too close to the 
field line or at the wrong angle relative to its equilibrium 
position. A careful launching of the signal can eliminate all 
these initial oscillations even though they may grow later when 
the ray crosses the equator or upon reflection. 

In Figures 10 and 11, the rectangular plots show one 
reflection each for the cases L = 1.63 and L = 1.53, respectively, 
while the corresponding polar plots show several round trips. For 
these two cases, the peak enhancement is 6 percent at the point of 
launch. 

Figure 12 shows how a small peak enhancement is sufficient 
(3.0 percent) to trap the rays if the signal is properly launched. 
The physical conditions are the same as for the case of Figure 9. 

Figure 13 relates to an L = 2 field line and shows the 
manner in which a ray typically escapes when the oscillations 
grow too large, even with sufficient peak enhancement. A fre- 
quency of 2.00 MHz is used for this case. 


National Aeronautics and Space Administration 
Electronics Research Center 

Cambridge, Massachusetts, July 1968 
188-39-01-01 
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Figure 10. - Ray tracing along the L = 1.63 field-line for 
6.0 percent peak enhancement at a geocentric distance of 
9231.57 km (frequency = 1 MHz, extraordinary wave) 
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Figure 11. - Ray tracing along the L = 1.53 field line for 
6.0 percent enhancement at a geocentric of 9455.50 km 
(frequency = 1 MHz, extraordinary wave) 
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Figure 12. - Ray tracing along the L = 1.76 field line for 
only 3.0 percent peak enhancement of electron density at 
a geocentric distance of 10,026.00 km (frequency = 1 MHz) 
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